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ABSTRACT 



Context. X-ray reflection off the accretion disc surrounding a black hole, together with the associated broad iron Ko- line, has been 
widely used to constrain the innermost accretion-flow geometry and black hole spin. Some recent measurements have revealed steep 
reflection emissivity profiles in a number of active galactic nuclei and X-ray binaries. 

Aims. We explore the physically motivated conditions that give rise to the observed steep disc-reflection emissivity profiles. 
Methods. We perform a set of simulations based on the configuration of a possible future high-resolution X-ray mission. 
Computations are carried out for typical X-ray bright Seyfert-1 galaxies. 

Results. We find that steep emissivity profiles with q ~ 4 - 5 (where the emissivity is e(r) oc r~ q ) are produced considering either 
i) a lamp-post scenario where a primary compact X-ray source is located close to the black hole, or ii) the radial dependence of the 
disc ionisation state. If both effects are taken into account, emissivity profiles as steep as q ~ 7 can be obtained from X-ray spectra 
modelled via conventional reflection models. We also highlight the role of the reflection angular emissivity: the radial emissivity index 
q is overestimated when the standard limb-darkening law is used to describe the data. 

Conclusions. Very steep emissivity profiles with q > 7 are naturally obtained by applying reflection models that take into account 
radial profile £(r) of the disc ionisation induced by a compact X-ray source located close to the central black hole. 

Key words. Black hole physics - Accretion, accretion discs - Relativistic processes - Galaxies: nuclei 
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1. Introduction 

The innermost black-hole accretion discs of active galactic nu- 
clei (AGNs) and black hole binaries may be revealed through 
their X-ra y radiation released as either the thermal radiation of 
the disc dZhang et all 1 19971: iMcClintock & Remillardll2006l) or 
the result of the inverse Compton scattering of the thermal pho- 
tons on the rel ativistically moving electrons in the hot coron a 
above the disc dShapiro et alj|1976t lHaardt & Maraschilll991l) . 
While the first is relevant only to stellar-mass black holes where 
the accretion disc is heated to very high temperatures (around 
10 7 K), the second is common for any accreting black hole over 
the entire possible mass scale. 

Some fraction of the photons scattered in the corona reflect 
off the disc surface before reaching the observer. This so-called 
reflection spectrum is characterised by a Compton hump at en- 
ergies of around 20 - 40 keV, fluorescent lines of which the iron 
Ka line at 6.4-6.97 keV is the most prominent, and a soft excess 
below 1 keV. The overall spectral shape depends on the ionisa- 
tion of the disc surface. The spectrum is then smeared by the 
relativistic effects including energy shift, aberration, and light- 
bending. The inner disc radius, inclination angle, average ioni- 
sation state, and reflection-emissivity radial profile can in princi- 
ple be determined from the resulting X-ray spectral shape. The 
inner d isc radius plays a particula rly important role. As shown 
e.g. by Reynolds & Fabian (2008), the inner-disc reflection ra- 
dius can be associated with the innermost stable circular orbit, 
which only depends on the black hole spin. Hence, X-ray reflec- 
tion spectra can be used to estimate the black hole spin in both 



AGN and black hole binaries (see e.g.lRevnolds & Nowakl2 003; 
Miller et al. 2009; Brenneman & R eynolds 2009, and references 
therein). 

The black hole spin measurements are thus influenced by 
the geometry of the disc-illuminating corona and the local prop- 
erties of the disc that affect the re-processing and re-emission 
of the incident photons. Current relativistic kernels that are ap- 
plied to reflection (and/or iron line) models to include the rela- 
tivistic effects on the spectral shape (lLaorl 1991; Dovciak etaD 
120041: iBeckwith & Pond 12004 iBrenneman & Reynolds! T2006) 
are based on a series of simplified assumptions and, in partic- 
ular, they assume a single/broken power-law form of the radial 
reflection emissivity and an angular emissivity law defined by 
a simple analytica l formula, m ost frequently employing a limb- 
darke ning profile jLa or 1991). 

In ISvoboda et alJ (|2009), we pointed out that the choice of 
a particular angular emissivity law can affect the derived rela- 
tivistic parameters, including the black hole spin. In the present 
work, we extend our analysis to the radial emissivity law and, 
in particular, we investigate the physical conditions under which 
steep radial emissivity profiles are produced by considering the 
effects of (i) primary X-ray source location, (ii) radial disc ioni- 
sation profile, and (iii) angular emissivity law. 

The intrinsic disc emissivity is naturally expected to decrease 
with increasing distance, i.e. the reflection emissivity is 



e(r) = r 



(1) 



where q is the emissivity index that can be constant over all 

radii or a varying quantity. The thermal dissipation of the disc 

Send offprint requests to: J. Svoboda, email: jsvoboda@sciops.esa.int decreases as r~ 3 (Shakura & Sunv aevlll973HNovikov & Thornd 
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Il973l) . Therefore, the simplest assumption is postulating the 
same dependence for the reflection. The more energetic photons 
are injected into the innermost regions, and so, more intense ir- 
radiation of the disc occurs there. In addition, assuming a point- 
like X-ray source at height h on the disc axis, the irradiation of 
the disc in the absence of any relativisti c effect is proportional 
to (r 2 + h 2 y 3 ^ 2 oc r~ 3 , as shown e.g. bv lRevnolds & Begelman 
(1997). An emissivity profile with q — 3 is therefore consid- 
ered as standard, while steeper/flatter indices may need to be 
explained. 

Steep emissivity profiles have been reported from the anal- 
ysis of_fhe_2^aysp_ectra of other AGN, such as MCG -6-30- 
15 dFabian et al.ll2002t IVaughan & Fabianl 12004 IMiniutti et al l 
1 20071). 1H0 707-4 95 dFabian et al.1 120091; [Zoghbi et all 1201 Ot 
Wilki ns & Fabian|l201 \i iDauser et aLJl2012J) . and IRAS 1 3224- 
3809 dPonti et al. 2010), as well as b lack hole bi naries, such a s 
XTE J 1650-5 00 (iMiniutti eTaill2004l) . GX 339-4 (lMilleril2007l) . 
and CygX-1 dFabian et al.l 1201 21) . The measured indices reach 
values up to q « 7. 

To provide a p hysical pictu r e of th e steep radial emissivity 
in MCG -6-30-15. IWilms et al.1 (l200lh invoke strong magnetic 
stresses that should act in the innermost region of the system. 
This should correspond to the enhanced dissipation of a consid- 
erable amount of energy in the accretion disc at small radii. If 
the magnetic field lines thread the black hole horizon, the dissi- 
pation could be triggered by the magnetic extraction of the black 
hole rotational energy, perh aps via the Blandford-Znajek effect 
(Blandford & Znaiek 1977), but it could be also supplemented 
by a rather efficient slowin g of the rotation, as also seen in recent 
GRMHD simulations (e.g. lPenna et al.ll2010h . The efficiency of 
the competing proces s es still needs to be assessed. 

iMartocchia et al.l (|T000) examined whether the required 
steep emissivity law as well as the predicted equivalent width of 
the cold reflection line of iron and the Compton reflection com- 
ponent can be reproduced in a phenomenological (lamp-post) 
model where the X-ray illuminating source is located on the 
common symmetry axis of the black hole and the equatorial ac- 
cretion disc. They suggested that the radial emissivity function 
of the reflection component steepens when the height parame- 
ter of the primary irradiation source decreases. The enhanced 
anisotropy of the primary X-rays was identified as a likely agent 
acting in this process. The emissivity in the XMM-Newton spec- 
trum of MCG -6-30-15 w as successfully reproduced by adopting 
this lamp-post geometry (jMartocchia et al.ll2002t IMiniutti et al.l 
l2003HNiedzwiecki & Zvckill2008l) . 

To explain the steep radial emissivity, we explore several 
simple test models and analyse them with the simulated data. 
The paper is organised as follows. Sect. [2] describes the data 
generation. The effects of the lamp-post geometry of the corona, 
the angular directionality, and the radially stratified ionisation on 
the radial emissivity are discussed in Sects. [3] [4] and [5] respec- 
tively. An example of the combined effect is presented in Sect. [6] 
The obtained results are discussed in Sect. Qand summarised in 
Sect. E 

2. Data simulation 

Our goal is to build a series of physically motivated reflection 
models and to examine the radial emissivity law that would 
be inferred from deep, high-quality X-ray spectral observations 
of AGN and/or black hole X-ray binaries with either current 
or future X-ray facilities. We build theoretical reflection spec- 
tral models that are then put into X-ray spectral simulations 
to produce high quality X-ray data. The simulated data are 



then modelled with the currently available disc-reflection mod- 
els and standard spectral tools. To evaluate the possibility of 
measuring the simulated effects with future X-ray missions, we 
adopted a preliminary response matrix for the Athena X-ray 
mission (NandriHoIi]). This ensures that our results are not 
driven/limited by the signal-to-noise that can be reached with 
current X-ray observatories. 

We performed our spectral analysis within the 2-10keV en- 
ergy range where one of the most prominent reflection feature, 
the iron Ka line, occurs. The preliminary response matrixQ 
was re-binned by a factor of ten between channels 2700 and 
8800 (2-10keV). The other channels were not used. This pro- 
vided us with sufficient spectral resolu tion and incre ase in the 
computational speed. We used Xspec (lArnaudlll996l) . version 
12.6.0ab for the spec tral fitting. We used the most recent ver- 
sion of the KY code dDovciak et al.ll200l . which includes the 
lamp-post geometry (Dovciak et al., in prep.). The flux of the 
model was chosen to be similar t o that of bright Sey fert galax- 
ies, i.e. w 3 x 10 ergs cm - 2 s 1 dNandra et alJl2 007). The sim- 
ulated observation time was 100 ks, which is comparable to the 
average exposure time of AGNs observed with the current X-ray 
satellites. 

3. Lamp-post scheme 

We first investigate how the radial emissivity depends on the ge- 
ometry of the corona. If the corona is localised the illumina- 
tion of the disc decreases with the growing distance from the 
source in a way determined by the position of the corona and 
by the gravitational field of the central black hole. The configu- 
ration when the corona is very compact and located just above 
the black hole, known also as the lamp- post scheme, was stud- 
ied as a simple d i sc-cor ona scenario by iMatt et all d!991l) and 
iGeorge & Fabianldl99ll) . In a physical picture, the source above 
the black hole can be imagined, e.g., as the base of a jet. 

In this scenario, the irradiation far from the source radially 
decreases as r~ 3 . In the central region, the relativistic effects in- 
fluence the disc illuminati on, thus shape the reflec tion spectra of 
black hole accretion discs dMiniutti & Fabianll2004l) . As a result, 
the different parts of the disc are irradiated with different intensi- 
ties, making the emissivity profiles in reflection models distinct 
from the standard value of q = 3. If the height of the source is 
sufficiently close to the black hole event horizon, light bending 
implies higher irradiation of the innermost region compared to 
the outer parts of the disc. 

The exact profile of the radial emissivity depends on the ge- 
ometrical properties of the source. Different c ases of axial, orbit- 
ing, je t, and extended sources were studied by Wilkins & Fabian 
(2012). The steepest profiles were obtained for point- like 
sources at small heights along the vertical axis. In the most ex- 
treme scenario, when the source is moving towards the black 
hole, the Doppler boosting might increase this effect. However, 
there is no observational evidence for such an inflow of the 
matter perpendicular to the disc plane, while outflows in the 
form of a jet o r a disc wind are observed in many sources (e.g . 
Merloni et al] | 2003[ iBlustin et all 120051: iKrongold et al.l l2007t 
Tombesi et atlboiOl) . 

In the following, we consider the stationary lamp-post source 
and investigate the radial emissivity profile of the disc reflection 
radiation for different heights of the source. The physical set-up 
of the model is a combination of the general-relativistic lamp- 
post scheme for an X-ray illuminated accretion disc near a ro- 

1 ftp://ftp.rssd.esa.int/pub/athena/0905201 l_Responses 
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Fig. 1. Contour plots of the spin a and the radial emissivity parameter q. The data were generated with the lamp-post model with 
the height h = 1.5 r g . The default value of the spin was a = 0.94, which is indicated by a dashed line in the graph. Different 
prescriptions for the angular emissivity were used: Top left: Angular emissivity from numerical calculations. Top right: Limb 
brightening. Bottom left: Isotropic. Bottom right: Limb darkening. The x 1 values corresponding to the best fit (indicated by a 
small cross) and to the lcr, 2<x, and 3<x levels are shown in the text legend. 



Table 1. Inner radial emissivity index q m and the break radius inferred for a different height and directionality in the lamp-post 
model. 

a = 0.94 
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tating black hole (Marto cchia et al.ll2000h and a self-consistent 
Monte Carlo scheme for the X-ray r eprocessing within the disc 
environment dGoosmann et al. 2006). We used local (re-) emis- 
sion table s that were compu ted by the radiative transfer code 
NOAR dDumont et al.l HOOO, Sect. 5) for the case of "cold" 
reflection (i.e. for neutral or weakly ionised matter). Photo- 
absorption, Compton scattering, and the fluorescent emission of 
the iron K line are considered. A stratified plane-parallel atmo- 
sphere irradiated by a power law with the photon index T — 1 .9 
is assumed. A large number of primary photons were sampled in 



the 2 - 300 ke V energy range to ensure sufficiently high quality 
statistics. At all local emission angles, the Poissonian noise level 
is much smaller than any relevant spectral feature. The computa- 
tions were done for various incident local emission angles, both 
polar and azimuthal. 

In the present paper, we employ a new implementation 
of the lamp-post sce nario that is based on the KY package 
dDovciak et al.l [20041 a more detailed description of the new 
model is presented in the parallel paper Dovciak et al., in prep.). 
All the relativistic effects on the photon energies and trajecto- 
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Fig. 2. Contour plots of the radial emissivity and the break radius for the height h — 1.5 r g . The legend is the same as in Fig.Q] A 
horizontal dashed line representing the best-fit value from Table [T] q = 4.72, is in the graph only for comparison. 



ries, such as aberration and light bending, are incorporated in 
the model. Here, we study whether it is possible to approximate 
the radial emissivity of the lamp-post model by a simplified pro- 
file in the form of a broken power-law, as usually done in current 
modelling of the data. 

For this purpose, we generated data as described in Sect. [2] 
for a power-law model with photon index F = 1 .9 and the rela- 
tivistic line in the lamp-post geometry described above. The spin 
parameter was set to a — 0.94 (in the geometrised units where 
c - G — M — 1). The other parameters are the inclination angle 
of the disc i = 30 deg, which is a typical value for the inclination 
of Seyfert 1 galaxies dde La Calle Per ez et aTll2010h . the inner 
radius coinciding with the marginally stable orbit (r; n = r ms ) and 
the outer radius corresponding to r out = 400 r % , where r g = is 
the gravitational radius. We considered three cases of the source 
height, h = 1.5, 3, and 10 r g . 

We then fitted the radial emissivity profile resulting from the 
simulations with a broken power-law model in which the emis- 
sivity index is q m within a break radius r\> and <7 ollt = 3 at larger 
distances. As expected, we found that q m increases as the source 
height h decreases, while r\, decreases. Steep emissivity profiles 
in the inner regions (i.e. q\ B > 3) were only obtained for h < 3 r g , 
as summarised in Table[T]in the column labelled "numerical". 



4. Interplay between the radial and angular 
emissivity profiles 

When fitting the data, the local intensity of the re-processed ra- 
diation emitted from the disc is often assumed to be divided into 
two separate parts - the radial and angular dependence. The lat- 
ter one characterises the emission directionality. However, owing 
to the large rotational velocity of the disc and the strong gravity 
near the black hole, the photons that reach the observer are emit- 
ted under different angles at different locations. Therefore, the 
angular part of the emissivity depends on radius as well, and 
the above separation is invalid. The relativistic effects, namely 
aberration and light bending, cause the emission angle in the in- 
nermost region to always be very high (almost 90 degrees wit h 
respect to t he disc normal) - see A ppendix C in Dovciak (2004), 
or Fig. 3 in lSvoboda et al.l d2009l) . Although it is not an axisym- 
metric problem, the almost radial decrease in the emission angle 
is readily apparent, whi ch invokes the link between the radial 
and a ngular emissivity ( Bec kwith & Donell2004t [Svoboda et al.l 
I2009h . 

To show the impact of the directionality on the radial emis- 
sivity profile, we considered the cases of 

1 . Our numerical computations, 

2. Limb brightening /(// e ) « ln(l + fx' 1 ) dHaardtlll993h . 

3. Isotropic emissivity, 

4. Limb darkening 7(ju e ) * 1 + 2.06// e dLaorill99ll) . 
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where /u e is the cosine of the local emission angle. The simu- 
lated data were created with our numerical model of the direc- 
tionality integrated o ver all incident angles. The NOAR code 
(DumontetaU]2000) was employed in the calculations. Free- 
free absorption, the recombination continua of hydrogen- and 
helium- like ions, and both direct and inverse Compton scatter- 
ing were taken into account. 

Figures Q] and [2] show contour plots of the radial emissivity 
index and the spin, and the break radius, respectively, for dif- 
ferent prescriptions of the angular directionality. The lamp-post 
model with the source height at h = 1 .5 r g was used in the data 
simulation, while the fitting model employed a broken power- 
law emissivity. In the contour calculations, only the two param- 
eters of interest were allowed to vary. The others were fixed to 
their default or best-fit (in the case of break radius) values. The 
best-fit parameters are summarised in Table Q] 

The worst fit was obtained for the model with limb darken- 
ing for which significantly steeper values of the radial emissivity 
were required. This emissivity law is widely used in the reflec- 
tion models. However, it somewhat cont radicts several models 
of X-ray illuminated disc atmospheres dGhisellini et alj [ l 994 ; 
Zvcki & Czernvl Il994t iGoosmann et alJ |2006; Roz ahska et alj 
201 ll) . Its application causes an appreciable underestimation of 
the innermost flux. The fit therefore requires an artificially steep 
value of the radial emissivity index in order to compensate for 
this loss of counts from the central region where the emission 
angle is grazing. 



5. Radially structured ionisation of the disc 

The interplay between the radial and angular emissivity shows 
that the steep radial emissivity in the observational data might 
be caused by an invalid model assumption. Yet, there is another 
frequently used assumption in the reflection scenario that can 
contribute to this effect as well - the constant ionisation over the 
whole surface of the disc. The intensity of the disc irradiation, 
whether approximated by a (broken) power-law decrease or a 
lamp-post illumination in curved space-time, decreases with ra- 
dius. Therefore, the ionisation of the disc surface may respond 
accordingly (see e.g. lMatt et al.ll 1993h . 



The importance of the photoionisation of th e disc sur- 
face in AGNs was studied b y several authors dZvcki et alj 
[1991: INavakshin et al.1 120001 iBallantvne etail 1200 ll 120111) . 
iBallantvne et al.1 (120031) interpreted the X-ray reflection spec- 
trum of MCG -6-30-15 as being composed of a highly ionised 
reflection from the innermost region and a cold one from the 
outermost accretion disc. The radially dependent ionisati on in 
the existing data was also discussed by IZhou et a D j20Tlh . The 
photoionisation was also suggested as a possible explanation 
of the non-detection of the spectral imprin ts of the relativisti- 
cally smeared reflection in some sourc es (Reynolds et al. 2004; 
Svobo da et al]l2010HBhava"nT & Nandra ll201lt Brenneman et al. 
120121) . 

In general, the ionisation of the disc surface depends on 
several other physical quantities such as the density, vertical 
structure, and thermal heating (see e.g. INa vakshin & Kalhnanl 
12001b Rozahska et al 1 120021: Goosmann et al. 2007, and refer- 
ences therein). In particular, INavakshin et alJ d2000l) showed 
that considering hydrostatic and ionisation balance and radia- 
tive transfer in a plane-parallel geometry resulted in a cold core 
within the accretion disc and highly ionised outer layers. A de- 
tailed description of the disc ionisation is beyond the scope of 
this paper. Here, we simply assume that the radial dependence 
of the ionisation may be relevant, as a natural consequence of 
the radial dependence of the disc illumination by the primary 
radiation. 

If the photoionisation is the dominating factor determining 
the state o f the plasma the ionis ation parameter can be defined 
as (see e.g. lRoss & Fabianll 19931) 



47T.Fi, 



«H 



(2) 



where F mc is the incident (irradiation) flux and «h is the number 
density of hydrogen. Solar abundances are assumed for the other 
elements. Typically, X-ray spectra are fitted using only one re- 
flection component, i.e. assuming £ constant over the whole disc. 
However, in general, eq.|2]may be re-written as 



(3) 



with the radial dependence of the irradiating flux and the density. 

Figure [3] shows the radial dependence of the ionisation for 
different density profiles: constant disc density, Novikov-Thorne 
dNovikov & Thornd 1973b . and a profil e from general relativistic 
magneto-hydrodynamic simulations by Hawley & Krolik (2006) 
for isotropic irradiation decreasing like F\ nc oc r~ 3 . In all cases, 
the ionisation decreases with the radius according to a decreas- 
ing strength of the irradiation. Although an unknown density 
profile introduces an uncertainty in the radial ionisation descrip- 
tion, its effect is, however, not as strong a function of the radius 
as the irradiation profile itself. We therefore assumed a constant 
density over the whole disc surface in our more detailed analysis. 

To describe the ionisation profile, a specific ionisation pa- 
rameter has to be set at a given radius. In general, this normalisa- 
tion depends on the luminosity and the distance of the irradiating 
source, in addition to the density of the disc. It is therefore differ- 
ent for different kinds of sources, and moreover, depends on the 
variable accretion rate. The choice of the normalisation factor is 
thus somewhat arbitrary for our purposes. Here we assumed that 
the ionisation of the disc is £ = 10 ergs cm s~' at 30 r g since the 
neutral iron line is observed in most spectra of AGNs. 

We simulated a spectrum produced under these conditions by 
dividing the disc into discrete rings of fixed £ and summing the 
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Fig. 4. Plot of the 'complex' model. The top solid line shows the total model flux, while the dotted lines are its components. The sum 
of all reflection components is shown by the green solid line. The red dashed line represents the best-fit reflection model assuming 
single ionisation. See the main text for more details. 



~ 10 

i> 

CD 

V co 10" 1 
p. 

X p 

i 10" 2 
1.2 



o 
'a 



x_ 



1 



0.8 




700 



680 



660 



640 









/ 3a 





10 



3.5 



4.5 



5.5 



Energy [keV] 



Fig. 5. Left: Data created by the "complex" model (see Figure [4]i and their ratios with respect to the single-ionisation reflection 
model. Right: Dependence of the fit-goodness on the radial emissivity parameter of the single reflection model. The dashed line 
corresponds to a 3<x confidence level. 



spectra produced in each ri ng. To generate these spectra, we used 
a grid of reflionx models (Ross & Fabi anll2005h with a difFerent 
ionisa tion state convolved with the KY model dDovciak et alj 
2004), corresponding to different emission regions across the 
disc. For each ring, the generated reflection component has equal 
integrated flux and the radial emissivity is governed by the rela- 
tivistic convolution models. The innermost radius coincides with 
the marginally stable orbit, and the outer radius was set to 400 r g . 
The spin value was chosen to be a = 0.94, i.e. r ms ~ 2 r g . The in- 
clination angle was chosen to be 30 deg. The photon index of the 



primary power-law radiation and its normalisation were set to 
F = 2.0 and Kr = 10~ 3 . The standard value, q = 3, was adopted 
for the radial emissivity index. We used solar iron abundances. 

The total reflection fraction to the primary radiation corre- 
sponds to the case of an isotropically irr adiated disc, i.e. similar 
to tha t of R a 1 in the PEXRAV model (Magdziarz & Zdziarskil 
Such a constructed model, which we henceforth refer to 
as a "complex" model, is shown with its components in Figure|4] 
The straight dotted line in the plot corresponds to the power-law 
continuum irradiating the disc. The iron line around 6 keV may 
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Fig. 6. Contour plots of the radial emissivity parameter q and the break radius r\, (left) and the spin a (right), respectively. Data were 
created by a "complex" model with radially stratified ionisation, isotropic irradiation, and angular emissivity. A single-ionisation 
model with power-law radial emissivity was used to fit the data. The individual contours correspond to l<x, 2cr, and 3<x level, 
respectively. 



Table 2. Resulting parameter values of the single reflection 
model applied to the data simulated by a "complex" reflection 
model. 



parameter 


isotropic 


h = 1.5 r g 


spin 


0.86!™ 


0.94 + 0.02 


inner radial emissivity 


4 <:+0.6 


6 7 +0 - 9 

u - '-0.8 


ionisation [ergs cm s~'] 


250 +3 ° 


230^ 


fit goodness x 1 l v 


639/606 


623/606 



Notes. Two cases of assumed irradiation profile are considered, 
isotropic with the radial emissivity index q = 3, and a point source 
at the height h = 1.5r g . Default value of the spin was a = 0.94. 



be used as a diagnostic of the individual reflection components. 
The cold reflection from the outer disc (13.9 - 400 r g ) corre- 
sponds to the most prominent peak. The second largest peak is 
produced by the ring occurring at 6.3 - 13.9 r g . The components 
from the more central regions have more smeared profile and 
more red-shifted iron lines. The disc separation into individual 
rings was done according to the gradient of the ionisation pa- 
rameter, so the innermost rings are much narrower than the outer 
parts and thus contribute less to the total flux (the bottom lines in 
the graph). The total reflection component is shown by a green 
line in the plot. 

We generated the data in the same way as described in 
Sect. [2] We then fitted the data in the 2-10keV energy range with 
a model consisting of only a power-law continuum and a single 
reflection component with a broken power-law radial emissivity. 
The inner radial emissivity index, the break radius, the ionisa- 
tion, and the normalisation of the reflection model were the only 
parameters that were allowed to vary during the fitting proce- 
dure. 

The resulting single-ionisation model represents a relatively 
good fit to the simulated data with^ 2 /v = 639/606 ~ 1.05, but 
with significant discrepancies in the iron line and its edge. The 



model (indicated by a red line) is compared with the "complex" 
reflection in Figure [4] An apparent difference occurs at the iron 
line, where a peak from distant cold reflection is missing in the 
single-ionisation model. The data residuals from this model are 
shown in the left panel of Fig. [5] The fit improved by A^ 2 = 16 
when a Gaussian line was added to the model. The energy and 
the width of the additional line were found to be E = 6.45+q 
and cr = 0.13+°<g. The equivalent width is 41 + 17 eV. 

The best-fit values of the reflection model parameters are 
summarised in the middle column of Table [2] The value of the 
ionisation parameter is around 250 ergs cm s _1 . The radial emis- 
sivity index is required to be significantly steeper, q > 4, as it 
is also clearly visible in the right panel of Fig. [5] Figure [6] (left 
panel) shows a contour plot between the radial emissivity param- 
eter and the break radius. A preferred value of the break radius 
is larger than 10r g , which demonstrates the significance of the 
emissivity steepening (i.e. that it is not related only to a nar- 
row innermost ring). The fit is insensitive to larger values than 
about 20r g . A contour plot between the radial emissivity parame- 
ter and the spin is shown in the right panel of Figure[6] A mutual 
degeneration of these two parameters is evident from the plot. 
However, significantly steeper values of the radial emissivity are 
required for any reasonable value of the spin. 



6. Steep radial emissivity due to the combined 
effect 

We have presented three different effects that may account for 
a steep radial emissivity. However, they are not completely in- 
dependent and the most likely scenario is that they occur to- 
gether. We therefore simulated the data by considering all these 
effects together. We considered a lamp-post scenario with a very 
low height, h = 1.5r g . We calculated the ionisation profile of a 
constant-density disc produced by this irradiation taking into ac- 
count all the relativistic effects such as photon light-bending and 
aberration. For simplicity, we assumed an isotropic directional- 
ity, i.e. that the intensity of reflected radiation does not locally 
depends on the emission angle. 

The data were then fitted by a single-ionisation reflection 
model with the radial emissivity defined by a broken power-law 
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Fig. 7. As in Figure [5] but the data were created by a lamp-post model with the height h — 1.5 r g . In the right panel, we show a 
further overestimation of q by employing limb darkening instead of isotropic angular directionality used in the data simulation (red 
dashed line). 



with the index q as a free parameter within 30r g and q — 3 be- 
yond this radius. The spectrum i s reflection-dominated d ue to 
the strong light-bending effect dMiniutti & Fabian] 12004). The 
parameters of the primary power-law were fixed. The resulting 
model describes the data quite well with x 2 l v = 623/606 ~ 
1.03, again with some residuals in the iron-line energy band (see 
the left panel of Figure|7]). The best-fit values are summarised in 
the last column of Table [2] 

The radial emissivity is very steep, q « 7. We show the de- 
pendence of the fit goodness on this parameter in the right panel 
of Figure|7] where we compare two cases of assumed direction- 
ality. A value larger than 5 within the 3<x uncertainties is found 
even with the isotropic angular emissivity. The limb darkening is 
responsible for the further steepening. Figure[8]shows a contour 
plot between the radial emissivity and the spin. In contrast to the 
radial emissivity, the spin is found to be close to its default value. 

7. Discussion and conclusions 

We have addressed the steep radial emissivities that have been 
detected in the reflection components of the X-ray spectra of 
active galaxies and black-hole binaries. We investigated some 
possible explanations, and to this end performed several simu- 
lations to reveal the degeneracies of the radial emissivity with 
other parameters and the intrinsic assumptions of the relativis- 
tic reflection model. We realised that the steep radial emissivity 
may be explained by either (i) the geometrical properties of the 
disc-illuminating corona, (ii) the use of an improper model as- 
sumption about the angular directionality, or (iii) radially struc- 
tured ionisation. The first puts rather extreme requirements on 
the corona. It needs to be very bright and occur at a very low 
height above the black hole. The second is likely due to use of 
an improper prescription for the emission directionality in the 
black hole accretion disc. The last one is related to a probable 
radial dependence of the disc ionisation, which plays a signifi- 
cant role in the total shape of the reflection spectrum. 

7. 1 . Lamp-post scenario 

The steep radial emissivity may be related to the proper- 
ties of the disc-ill uminating corona as suggested before by 
IWilms et al.l d200ll) . The geometry of the emitting region cer- 



tainly plays a significant role. A very centrally localised source 
at a low height above the black hole horizon would irra- 
diate the disc mainly in its central region. The illumination 
in this area is greatly enhanced by t h e gravitational light- 
bending effect dMiniutti & Fabianl l2004t iNiedzwiecki & Zvckil 
I2008I; IWilkins & Fabiaij20JJl Dovciak et al., in prep.). 

To achieve steep radial emissivity, which is assumed to be 
proportional to the illumination, the source must be sufficiently 
close to the black hole. However, in this case the primary emis- 
sion has to be extremely bright because only a small fraction 
would overcome the strong gravitatio nal pull of the black hole 
and reach the observer (see Fig. 2 in Dovc iak et al.ll20lTl) . The 
importance of these effects declines sharply with height. At 
heights h > 3 r g , the radial emissivity profile is similar to the 
simple power-law with the standard value (q = 3). For even 
larger heights, the irradiation profile is more complicated (see 
Fig. 3 and 4 in iDov ciak et al. 201 TJ). It decreases steeply only 
very close to the black-hole horizon, then becomes rather flat 
(q < 3) when still in the inner parts of the disc and finally reaches 
the standard value far from the centre. 

7.2. Angular directionality 

For the angular emissivity, the limb darkening law is frequently 
used. Several simulations, however, suggest that the directional- 
ity is the opposite of limb darkening (see e.g. Rozanska et al.l 
2011, and references therein). The emission angle in the in- 
nermost region of the disc is always very high owing to strong 
aberration. The flux contribution from this region is therefore 
underestimated by models with limb darkening if the angular 
emissivity is indeed different. This effect could lead to an ap- 
proximately 20 % overestima tion of either the sp in or the inner 
radial emissivity parameter. lSvoboda et al. (2009) reanalysed the 
XMM-Newton observation of MCG -6-30-15 and showed that 
the radial emissivity might be a more sensitive parameter to the 
angular directionality than the spin. This is especially true when 
the spin value itself is very high (close to one). 

7.3. Radially structured ionisation 

We also discussed the impact of the probable radial dependence 
of the disc surface ionisation. The disc illumination by a corona 
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Fig. 8. Contour plot of the radial emissivity parameter q and the 
spin a. The data were created by a lamp-post model with the 
height h — 1 .5 r g , radially stratified ionisation, and isotropic an- 
gular emissivity. A single-ionisation model with limb darkening 
and power-law description for the radial emissivity was used to 
fit the data. The best-fit parameters are indicated by a small cross 
within the contours. 



is commonly assumed to be stronger in the innermost regions. 
We therefore assumed that the ionisation is higher in the inner- 
most region and decreases with radius. We did not consider other 
aspects that affect the ionisation structure of the disc, such as the 
density profile, vertical structure, and thermal processes (the last 
one being especially relevant for the stellar-mass black hole bi- 
naries). We simply assumed that the radial dependence of the 
irradiation is the dominating effect in determining the ionisation 
of the disc surface. This is especially true when the lamp-post 
scenario for a source of low height is considered. 

We note, however, that the density may still play an impor- 
tant role around the marginally stable orbit. While the Novikov- 
Thorn e profile of the density diverges there dNovikov & Thornel 
fl973h . i.e. the ionisation would decrease to zero, the accretion 
disc may extend below the marginally stable orbit owing to the 
presence of a magnetic field, as in the GRMHD simulations of 
Hawl ev & Kro lik (2006). The density does not increase to an in- 
finite value, bu t instead decreases close t o the marginally stable 
orbit (see also Rey nolds & F abian 2008). This region was, for 
clarity, not shown in Figure [3] but we plan to investigate it in 
more detail in a follow-up paper. 

By fitting the simulated data, we realised that the radial 
decrease in the disc ionisation may account for the radial- 
emissivity steepness equally as well as the assumption of the 
centrally localised corona. For the case of an isotropically illu- 
minated disc, we obtained radial emissivities of q as 4 - 5, i.e. 
values similar to those for the lamp-post irradiation of a cold 
disc with the height h = \.5r„. This is due to the different shapes 
of the reflection model for different ionisation parameters. The 
softer, more ionised reflection comes from the innermost part 
of the disc. When a simplified model with a single ionisation is 
used to fit the data, it may lead to a significant underestimation 
of the flux from the innermost regions, which is then artificially 
compensated for by a steep value of the radial emissivity profile 
in the model. 

The main difference between the single-ionisation model and 
the data created by a complex ionisation occurs at the iron-line 
energy band and its edge (see Figs. [5] and [7J . While the shape 



of the "complex" reflection continuum may be mimicked by 
a modest value of the ionisation and steep radial emissivity, the 
iron line peak due to reflection from the cold distant parts of 
the disc is not incorporated in the simplified model. An addi- 
tional emission line is then needed to model the residuals. In 
our test case, the equivalent width of such a line was found to 
be about 40 e V. Se veral AGN spectra, such as MCG -6-30-15 
dFabian et alJ l2002). have been found to contain, together with 
the broad disc component, a narrow iron line associated with 
reflection from a distant torus. A typical equivalent width is 
~ 100 eV. The reflection from the outer parts of the accretion 
disc might provide a significant contribution to this line when 
the single-ionisation reflection model is applied to the disc-line 
modelling. 

In our simulations, we considered the radial dependence 
of the ionisation. However, a similar effect on the emissivity 
profile may also be caused by a vertically structured accretion 
disc where the core is cold but the o u ter layers are hot, i.e. 
strong ly ionised (Navak shin et al.ll2000h . lNavakshin & Kal lman 
(2001) showed that in the case of AGNs, when the illuminat- 
ing flux is much higher than the thermal radiation of the disc, 
the hot skin of the disc is completely ionise d and most of the re- 
processing occurs in the disc core. However. lDauser et al .l (f2012h 
revealed the significant presence of highly ionised reflection to- 
gether with relatively cold reflection from the disc in the X-ray 
spectrum of a narrow-line Seyfert 1 galaxy 1H0707-495. The 
double-reflection model is shown in their Figure 6. Owing to 
their different shapes, the ionised reflection contributes more to 
the red wing of the observed iron line than the cold component. 
Neglecting the ionised component would therefore result in the 
measurement of a steeper radial emissivity index. 

Finally, we note that the simulations with the radially chang- 
ing ionisation for the isotropic irradiation were done with the as- 
sumption of a "standard" value for the radial emissivity, q — 3. 
However, non-thermal coronal emission does not necessarily 
need to behave in the same way as the thermal dissipation of 
the disc. The interaction between the disc and the corona is 
more complica ted, including the radiation and magnetic pro- 
cesses (see e.g. | Haa rdt & Maraschi 1991; Czernv & Goosmannl 
2004; IGoosmann et alj |2006; Rozarisk a et al.ll201ll) . In particu- 
lar, when the magnetic field is considered, the intrins ic profile 
might already be as steep as r 5 dKawanaka et al.ll2.005h . Further 
investigation of these poorly understood disc-corona interactions 
is therefore desirable. 

7.4. The combined effect 

All the effects that we have discussed are not independent of 
each other. Irradiati on by a compact centra lly localised source is 
very anisotropic dMiniutti & F abian 2004). The central regions 
are more illuminated and thus, more likely to be ionised than 
the outer parts of the accretion disc. The effect of a hypotheti- 
cal radial stratification of the disc ionisation is stronger than in 
the case of isotropic irradiation. We therefore performed a final 
simulation including both effects together. We simulated the data 
with the lamp-post geometry and the ionisation stratification cal- 
culated from the theoretical irradiation profile (assuming con- 
stant density). We fitted the generated spectrum with a single- 
ionisation model with a broken power-law for radial emissivity. 
As expected, we obtained significantly steeper radial emissivity, 
q as 7 (see the right panel of Tableland Figure[8]). 

Angular directionality plays an important role in the to- 
tal emissivity. We showed that the model with limb darken- 
ing overestimates the index of the radial emissivity profile by 
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5 - 10% (see the right panel of Figure |7). The effect might be 
even stronger because our numerical calculations of local emis- 
sivity suggest that there is a limb-brightening effect, while we 
conservatively used isotropic directionality in our simulations. 
Moreover, in the case of anisotropic illumination by a lamp-post 
source, the emissivity also depends on the incident angle with 
a similar limb-brightening profile (Dovciak et al., in prep.). The 
highest output of the reflected emission is thus obtained when 
both angles are grazing (perpendicular to the disc normal). This 
happens especially very close to the black hole owing to the large 
aberration that is caused by the extreme velocities of the mat- 
ter in the inner disc. Hence, the flux from the innermost region 
might be significantly enhanced and thus the radial-emissivity 
index overestimated. 

8. Summary 

Very steep radial emissivity profiles with q ~ 7 have been previ- 
ously claimed from the analysis of the disc reflection component 
in the X-ray spectra of some AGNs and black hole binaries. We 
have shown here that these profiles with q ~ 4-5 can be a natural 
consequence of either i) a compact centrally concentrated X-ray 
corona located 1 — 2 r g away from the central black hole, or ii) 
the radial dependence of the disc ionisation state. By combining 
the two effects, emissivity indices as steep as q ~ 7 have been 
obtained, which strongly suggests that objects such as 1H 0707- 
495 and IRAS 13224-3809 are indeed dominated by a centrally 
concentrated X-ray corona whose radially-dependent disc irra- 
diation induces a ionisation profile £(r) + constant at the disc 
surface. We have also demonstrated that an additional steepen- 
ing can be simply an artifact of improper usage of the emission 
directionality in the relativistic models. 
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